Screening of a Moving Parton in the Quark-Gluon Plasma 
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The screening potential of a parton moving through a 
quark-gluon plasma is calculated using the semi-classical 
transport theory. An anisotropic potential showing a min- 
imum in the direction of the parton velocity is found. As 
consequences possible new bound states in the quark-gluon 
plasma and J/ij) dissociation are discussed. 

Screening of charges in a plasma is one of the most 
important collective effects in plasma physics. In the 
classical limit in an isotropic and homogeneous plasma 
the screening potential of a point-like test charge Q at 
rest can be derived from the linearized Poisson equation, 
resulting in Debye screening. In this way the Coulomb 
potential of a charge in the plasma is modified into a 
Debye-Hiickel or Yukawa potential [1] 
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with the Debye mass (inverse screening length) mr> (fr = 
c = ks = 1). 

A special kind of plasma is the so-called quark-gluon 
plasma (QGP), where the electric charges in a plasma 
are replaced by the color charges of quarks and gluons, 
mediating the strong interactions between them. Such 
a state of matter is expected to exist at extreme tem- 
peratures, above 150 MeV, or densities, above about 10 
times nuclear density. These conditions could be fulfilled 
in the early Universe for the first few microseconds or in 
the interior of neutron stars. In accelerator experiments 
high-energy nucleus-nucleus collisions are used to search 
for the QGP. In these collisions a hot and dense fire- 
ball is created which might consist of a QGP in an early 
stage (less than about 10 fm/c) [2]. Since the masses of 
the lightest quarks and of the actually massless gluons 
are much less than the temperature of the system, the 
QGP is an ultrarelativistic plasma. To achieve a theo- 
retical understanding of the QGP, methods from quan- 
tum field theory (QCD) at finite temperature are adopted 
[3]. Perturbative QCD should work at high temperatures 
far above the phase transition where the interaction be- 
tween the quarks and gluons becomes weak due to a spe- 
cific property of QCD called asymptotic freedom. An 
important quantity which can be derived in this way is 
the polarization tensor describing the behavior of inter- 
acting gluons in the QGP. From the polarization tensor 
important properties of the QGP, such as the dispersion 
relation and damping of the plasma modes or the Debye 
screening of color charges in the QGP can be derived [4]. 



In the QGP the Debye mass of a chromoelectric charge 
follows from the static limit of the longitudinal polariza- 
tion tensor in the high-temperature limit [4], 
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where g is the strong coupling constant and n f the num- 
ber of light quark flavors in the QGP with m q <C T. 

The high-temperature limit of the polarization tensor 
corresponds to the classical approximation. For instance, 
it is closely related to the dielectric function following 
from the semi-classical Vlasov equation describing a colli- 
sionless plasma. E.g., the longitudinal dielectric function 
following from the Vlasov equation is given by [5-7] 
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(The only non-classical inputs here are Fermi and Bose 
distributions instead of the Boltzmann distribution.) 
Quantum effects in the Debye mass have been consid- 
ered using the hard-thermal-loop resummation scheme 
[8], dimensional reduction [9], and QCD lattice simula- 
tions [10]. The Debye mass and polarization tensor have 
also been computed in the case of an anisotropic QGP 
[11-13]. 

The modification of the confinement potential below 
the critical temperature into a Yukawa potential above 
the critical temperature might have important conse- 
quences for the discovery of the QGP in relativistic 
heavy-ion collisions. Bound states of heavy quarks, in 
particular the J/tjj meson, which are produced in the 
initial hard scattering processes of the collision, will be 
dissociated in the QGP due to screening of the quark 
potential and break-up by energetic gluons [14]. Hence 
the suppression of J/tp mesons have been proposed as 
one of the most promising signatures for the QGP for- 
mation [15]. Indeed, a suppression of J/ip mesons has 
been observed experimentally [16] and interpreted as a 
strong indication for the QGP formation in relativistic 
heavy- ion collisions [17]. 

In most calculations of the screening potential in the 
QGP so far, it was assumed that the test charge is at rest. 
However, quarks and gluons coming from initial hard pro- 
cesses receive a transverse momentum which causes them 
to propagate through the QGP [18]. In addition, hydro- 
dynamical models predict a radial outward flow in the 
fireball [19]. Hence, it is of great interest to estimate the 
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screening potential of a parton moving relatively to the 
QGP. Chu and Matsui [20] have used the Vlasov equa- 
tion to investigate dynamic Debye screening for a heavy 
quark-antiquark pair traversing a quark-gluon plasma. 
They found that the screening potential becomes strongly 
anisotropic. 

In the case of a non-relativistic plasma the screen- 
ing potential of a moving charge Q with velocity v fol- 
lows from the linearized Vlasov and Poisson equations as 
[21,22] 

, i» Q f ,o, exp [— ik ■ (f — vt)] 
4>{r,t;v) = / d 3 k — ^ 4 

It is easy to show that this expression reduces to the 
Yukawa potential in the case of small velocities, v = 
\v\ <C v t h, where Vth is the thermal velocity of the plasma 
particles. In the opposite case, v ^> v t h, the Coulomb po- 
tential is recovered since a screening charge cloud cannot 
be formed for fast particles. 

The above equation (4) also holds in the case of a rel- 
ativistic plasma. We only have to use the relativistic ex- 
pression (3) for the longitudinal dielectric function. For 
small velocities, v — > 0, i.e. u) <C k, we obtain 
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from which again the (shifted) Yukawa potential results 
[22] 

4>(f,t;v) = ^ exp(-m D \f-vt\). (6) 
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It should be noted that the opposite limit v ^ v t h, lead- 
ing to a Coulomb potential in the non-relativistic case, 
cannot be realized in an ultrarelativistic plasma because 
the thermal velocity of the plasma particles is given by 
the speed of light, vth = c = 1. 

In the general case, for parton velocities v between 
and 1 , we have to solve (4) together with (3) numerically. 
Since the potential is not isotropic anymore due to the 
velocity vector v, we will restrict ourselves only to two 
cases, r parallel to v and r perpendicular to v, i.e., for 
illustration we consider the screening potential only in 
the direction of the moving parton or perpendicular to 
it. 

In Fig.l the screening potential <j>/Q in w-direction is 
shown as a function of r' = r — vt, where r = |r|, between 
and 6 fm for various velocities. For illustration we have 
chosen a strong fine structure constant a s = g 2 /(4ir) = 
0.3, a temperature T = 0.25 GeV, and the number of 
quark flavors nf = 2. The shifted potentials depend 
only on v and not on t as it should be the case in a 
homogeneous and isotropic plasma. For r' < 1 fm one 
observes that the fall-off of the potential is stronger than 
for a parton at rest. The reason for this behavior is the 
fact that there is a stronger screening in the direction of 



the moving parton due to an enhancement of the particle 
density in the rest frame of the moving parton. 

In addition, a minimum in the screening potential at 
r' > 1 fm shows up. For example, for v = 0.8 this mini- 
mum is at about 1 .5 fm with a depth of about 8 MeV. The 
occurrence of a minimum in the potential in the direction 
of the velocity is also observed in so-called complex plas- 
mas. Complex plasmas are classical, low-temperature 
plasmas containing particles with a diameter of a few 
microns [23]. These particles are charged in the plasma 
by collecting electrons on their surface. In the presence 
of an ion flow the positively charged ions are deflected by 
the microparticles, leading to an anisotropic distortion of 
the Debye sphere by an enhancement of the ion density 
in front of the microparticlc. This positive charge cloud 
leads to an attraction between microparticles in the di- 
rection of the ion flow and the formation of string like 
structures [24], observed in experiments. 

A minimum in the screening potential is also known 
from non-relativistic, complex plasmas, where an attrac- 
tive potential even between equal charges can be found 
if the finite extension of the charges is considered [25]. A 
similar screening potential was found for a color charge at 
rest in Ref. [26], where a polarization tensor beyond the 
high-temperature limit was used. However, this approach 
has its limitation as a gauge dependent and incomplete 
(within the order of the coupling constant) approxima- 
tion for the polarization tensor was used [27] . Obviously, 
a minimum in the intcrparticle potential in a relativistic 
or non-relativistic plasma is a general feature if one goes 
beyond the Debye-Hiickel approximation by either tak- 
ing quantum effects, finite velocities, or finite sizes of the 
particles into account. 

Note that Chu and Matsui [20] did not report the exis- 
tence of a minimum in the potential of a quark traversing 
the QGP. However, in their Fig. 1(d) a negative value of 
the potential of a fast quark (v — 0.9) in the direction 
of the quark velocity is shown. Since the potential has 
to tend to zero for large distances, this implies a mini- 
mum in the screening potential. This minimum was not 
found because the potential was plotted only for the lim- 
ited range < r < l/ra<j — 0.35 fm for our choice of the 
parameters. 

The minimum could give rise to bound states, e.g., of 
diquarks, if thermal fluctuations do not destroy them. 
The two-body potential, associated with the dipole fields 
created by two test charges Qi and Qi at f\ and f*2 with 
velocities v\ and V2, can be written as 

- r 2 ,vi - V2,t) = 
Q1Q2 f J exp [ik - (-(ri - r 2 ) - {vi -#2)*)] 
4tt 2 J j k 2 Re[ei{cj = k-v 1 ,k)} 

exp [ik ■ ((n - r 2 ) - (jft - v 2 )t)] \ ^ 
k 2 Re[ei(uj = k ■ V2,k)] I 

For comoving quarks, V\ — V2, this two-body potential 
reduces to the one-body potential (4) , showing the attrac- 



2 



tion between the quarks which could give rise to a bound 
state. Colored bound states, e.g., diquarks, of partons 
at rest have also been claimed by analyzing lattice data 
[28]. For a quark-antiquark system, where Q1Q2 < 0, 
the two-body potential is inverted, showing a maximum. 
This may lead to short living mesonic resonances and an 
enhancement of the attraction between quarks and anti- 
quarks of mesonic states moving through the QGP. 
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Fig. 1 : Screening potential parallel to the velocity of the 
moving parton in a QGP as a function of r — vt (0 to 6 
fm) for v = 0, 0.5, and 0.8. 

The details of the potential, e.g., the depth of the min- 
imum, depend on the choice of the parameters, such as 
the coupling constant. For a value a s ~ 0.3, as it is 
typical for the temperature reachable in heavy-ion col- 
lisions, the semi-classical approach corresponding to the 
weak coupling limit might not be reliable. Quantum ef- 
fects and collisions between plasma particles are impor- 
tant at those temperatures and will change the dielectric 
functions and dispersion relations. Within the transport 
theoretical approach collisions can be considered, for ex- 
ample by using the relaxation time approximation [29]. 

Non-abelian effects (beyond color factors, e.g., in the 
Debye mass) will be important at realistic temperatures. 
Unfortunately they cannot be treated by the methods 
used here and are therefore beyond the scope of this work. 
However, as we discussed above, also in a complex plasma 
which is a strongly coupled plasma as it is also the case 
for the QGP close to the critical temperature, there is 
an attraction between the microparticles in the presence 
of an ion flow. This appears to be a general feature of 
weakly as well as strongly coupled plasmas. Therefore we 



do not expect a qualitative change of the screening po- 
tential due to non-pcrturbative and non-abelian effects. 

The results for f perpendicular to v are shown in 
Fig. 2, where the potential is shown function of 

\r — vt\ = \J r 2 + v 2 t 2 between 0.1 and 1 fm. Here we con- 
sider only the case t = since at t > and v > there is 
no singularity in the potential due to y/r 2 + v 2 t 2 > for 
all r. Hence the potential is cut-off artificially at small 
distances if plotted as a function of \/r 2 + v 2 t 2 . In con- 
trast to the parallel case (Fig.l) the fall-off of the poten- 
tial at larger values of v is less steep, i.e. the screening 
is reduced as it is expected since the formation of the 
screening cloud is suppressed at large velocities. Also no 
minimum in the potential is found. 
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Fig. 2: Screening potential perpendicular to the velocity 
of the moving parton in a QGP as a function of \f — vt\ 
(0.1 to 1 fm) for v = 0, 0.5, and 0.8. 

Summarizing, we have calculated the screening poten- 
tial of a color charge moving through the QGP from semi- 
classical transport theory corresponding to the high- 
temperature limit. As in Ref. [20] we obtained a strongly 
anisotropic screening potential. The screening is reduced 
in the direction perpendicular of the moving parton but 
increased in the direction of the moving parton, which 
may lead to a modification of the J/ip suppression. In 
addition, we found a new feature of the screening poten- 
tial of a fast parton in a QGP: a minimum in the poten- 
tial shows up which could give rise to bound states of, 
for example, diquarks if not destroyed by thermal fluctu- 
ations. For a quark-antiquark pair this minimum turns 
into a maximum which could cause short living mesonic 
resonances. Combining the effect of reduced screening 
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in perpendicular direction and the presence of a maxi- 
mum in parallel direction we expect a stronger binding 
of moving J/ip mesons than of J/ip mesons at rest with 
respect to the QGP. The consequences, e.g, for the J/ip 
yield should be investigated in more detail using hydro- 
dynamical models or event generators for the space-time 
evolution of the fireball. Finally, let us note that our re- 
sults also applies to other ultrarelativistic plasmas such 
as an electron-positron plasma in Supernova explosions. 
In this case one simply has to replace the Debye mass 
ui d by eT, where e is the electron charge. 
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